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Abstract 

An elementary, albeit higher dimensional, argument is used to compute the area under the 
power function curve between and 1. 

1 Introduction 

The purpose of this paper is to give a geometric proof that Jq x" dx = ^^q:^ n > 0. The 
proof, although quite simple, depends deeply on the geometry of higher dimensional space. The 
argument is designed to inform the reader, presumably an advanced undergraduate, about some of 
the geometrical configurations that can be found in higher dimensional space. Thus, our point of 
view is to take a familiar concept — that of area under a curve — and reveal it in a different light. 

In order to analyze higher dimensional geometry, we need to use rigid transformations of these 
spaces. Such transformations are found among the set of square matrices. Thus our proof uses 
basic linear algebra to check all the details. It is our hope that by thinking geometrically about 
the matrix manipulations, the reader can gain geometric intuition about higher dimensional space. 

Let us sketch the proof. Let x E [0, 1], then the quantity y = x" is the n-dimensional volume 
of an n-cube each of whose edge lengths is x. Such an n-cube can be found embedded in the 
unit (n + l)-cube, at height x with its edges parallel to the coordinate vectors ei = (1,0, ... ,0), 
62 = (0, 1, 0, . . . , 0), through e„ = (0, 0, ... , 1, 0). That is to say that for x between and 1, the 
cube rests in the hyperplane Xn+i = x in a standard position. As x varies from to 1, the union 
of these n-cubes forms a pyramid with n-cubical base. The area under the curve x" coincides with 
the (n + l)-dimensional volume of the pyramidal cone on the unit base. Our main task then will be 
to show that (n + l)-of these pyramids fill the unit (n + l)-cube. To see this, we will show that the 
pyramid intersects a simplex dual to the diagonal line from (0, 0, ... , 0, 0) to (1,1,..., 1, 1) in the 
barycenters of some of its lower dimensional faces. The simplex has a cyclic symmetry of degree 
(n + l)-and the rotations can be applied to the pyramid so that its orbits fill the (n + l)-cube. 
Specifically, there is a symmetry of the unit (n+ l)-cube of order n + 1 which gives an action of the 
cyclic group, on the unit (n + l)-cube. The orbit of the pyramid under this action consists of 

(n+1) copies of the pyramid, which tesellate the (n + l)-cube. Thus the pyramid takes up l/(n + 1) 
of the volume thereof. Since the volume of the pyramid coincides (up to choice of unit) with the 
integral, x" dx, this gives the result. 
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In the preceding paragraph, a lot of terminology was introduced without explanation. Before 
we embark on the full higher dimensional journey, we first run through the idea in dimension 3. 
Then we examine the ideas in 2-dimensions. After those examples are fully worked, we introduce 
the general terminology in higher dimensional space, and describe the intersections between certain 
subspaces in higher dimensional space. Then we describe how the pyramidal shape can be rotated 
to fill the cube. We discuss filling hyper-cubes with simplices, we give some illustrations, and we 
give the Mathematica code that can be used to generate pictures of the pyramidal sets. 
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2 Three dimensions 

Here is the second most easy case. The quantity y = represents the area of a square all of whose 
edge lengths are x. We arrange a set of these squares in a unit cube as follows. The unit square 
lies on the top of a 1 x 1 x 1 cube as the set Si = {(u, v, 1) : < u < 1 & < w < 1}. The square 
of side length x, and area is found as Sx = {{u,v,x) : < u < x & 0<f< x}. The union 
of these squares Uo<i<i5'x forms a pyramid with a square base, height 1 and volume x^dx. The 
pyramid has two faces that are right isosceles triangles with sides of length 1, 1 and ^/2. For either 
of the remaining two right triangular faces, the hypotenuse is of length ^/3, a side of length 1, and 
one of length ^/2. The pyramid that we are considering, then is the intersection of all the convex 
sets in 3-space that contain the vertices (0,0,0), (1,0,1), (0,1,1), (1,1,1), and (0,0,1). 

We call the edge from (0, 0, 0) to (1, 1, 1) the main diagonal. It can be given in parametric form 
as i{t) = {t,t,t) for t € [0, 1]. The main diagonal intersects an equilateral triangle in the cube at 
the centroid of the triangle. 

The equilateral triangle in question has as its vertices (1,0,0), (0,1,0), and (0,0,1). Observe 
that the intersection of the triangle with the diagonals of the pyramid of length \/2 occur at 
the points (1/2,0,1/2), and (0,1/2,1/2) while the big diagonal intersects at (1/3,1/3,1/3). The 
triangle can be given as the intersection of the cube with the plane x + y + z = 1. Each of 
the little diagonals can be parameterized, and these intersections can be determined. The entire 
pyramid intersect the triangle in a kite shape with vertices (0, 0, 1), (1/2, 0, 1/2), (0, 1/2, 1/2), and 
(1/3, 1/3, 1/3). For convenience, we indicate the kite as D2[3] and write 

D2[3] = H{{{0, 0, 1), (1/2, 0, 1/2), (0, 1/2, 1/2), (1/3, 1/3, 1/3)}); 

the subscript 2 indicates that the kite is two dimensional, and the [3] indicates that the third 
component of each vertex is non-zero. The letter H indicates the convex hull which will be defined 
in general in Section [l] 
Under the rotation 

(1,0,0)^(0,1,0)^(0,0,1)^(1,0,0), 
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Figure 1: Tiling the three dimensional cube with pyramids 



the kite L>2[3] rotates to I?2[l] = i^({(l, 0, 0), (1/2, 1/2, 0), (1/2, 0, 1/2), (1/3, 1/3, 1/3)}), and this 
in turn rotates to L>2[2] = -H'({(0, 1, 0), (0, 1/2, 1/2), (1/2, 1/2, 0), (1/3, 1/3, 1/3)}). These three kites 
tile the equilateral triangle. The restriction of the rotation to the pyramid gives two more pyramids. 
The three of these tile the cube. Under the rotation, the square base rotates first to the x = 1 face 
and then to the y = 1 face of the cube. The y = face of the cube is tiled by two isosceles right 
triangles as are the faces x = 0, and z = 1. The illustration in Fig. Ogives all the details. 

What we have shown then is that the integral Jq dx = 1/3. The quantity on the left side 
represents the volume of the pyramid, and the quantity on the right represents 1/3 of the cube. 

3 The two dimensional case 

The first sentence of the preceeding section was a teaser. If the case n = 2 is the second most easy 
case, what is the easiest case? For n = 1, we arrange horizontal line segments of length x at height 
y = x m the plane. The union of these segments is the isosceles right triangle with vertices (0,0), 
(0, 1), and (1, 1). The hypotenuse of this triangle is the line segment £{t) = {t,t) for t E [0, 1]. The 
triangle clearly occupies half of the square, {(x, y) S : < x < 1 & < y < 1}. Moreover, we 
can reach the other half of the triangle by reflecting in the line i{t). Fig [2 contains the details of 
the construction. The analogue of the triangle is the line segment on x + y = 1 that lies in the first 
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Figure 2: Tiling the square with triangles 

quadrant (so > x,y). 

We remark here on an interesting difference between dimensions 2 and 3. In dimension 2 the 
symmetry of the square that takes the lower triangle to the upper one is a reflection in the diagonal. 
In dimension 3 it is a rotation. If we extend the former situation to all of 3-dimensions, then we can 
think of the reflection in the line as having been achieved by a rotation in 3-space. In the general 
situation, the even dimensional cubes will be reflected. 

4 General notation 

4.1 The {n + l)-dimensional cube 

We begin by describing the unit (n+ l)-dimensional cube. This is given abstractly as the Cartesian 
product: 

[0,l]x[0,l]x...x[0,l] = [0,ir+^ 
(n+i) -factors 
or more concretely as a subset of R"+-'^ 

C"+i = {ixi,X2, . . . , Xn+i) G M"+^ : 0<Xj<l for all j = 1, 2, . . . , n + 1}. 

In R"+^ it is convenient to denote the standard basis vectors as follows: ei = (1,0, ...,0), 
62 = (0,1,..., 0), e„ = (0, 0, . . . , 1, 0), and e„_|_i = (0, 0, . . . , 0, 1). In general, throughout 
our discussion n is a fixed integer and n > 0. In the case that n = 0, the standard basis vector is 1. 
The important properties of a basis such as £^ = {ei, 62, ... , e„_|_i} are 

1. Any vector v £ R"+^ can be written as a linear combination — v = '}2^=i ^j^j where aj G M 
of elements in £; 

2. Such an expression is unique. That is, if v = Yl^=i ^j^j ~ X]j=i (^j^j^ then aj = /3j for all 
j = 1,2,..., n + 1. 

A set, S, is convex if any two points that are in the set can be joined by a line segment 
that is contained in S. If S" C M"^^, and P,Q £ S, then the line segment from P to Q can 
be given in parametric form as CpQ^t) = P + t{Q — P), for t £ [0,1]. Here we are using the 
vector space structure of M"+^ so that if P = (xi,X2, ■ ■ ■ ,Xn+i) and Q = {yi,y2, ■ ■ ■ lUn+i), then 
Q — P = {yi — xi,y2 - X2,.. ■ ,2/n+i - Xn+i)- A point on the line segment {ipqit) : t G [0, 1]} is 
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determined by vector addition as well. In the jth. coordinate, we have Xj+t{yj — Xj) = {l—t)xj+tyj. 
At t = 0, the point ^pq(O) = P and £pq(1) = Q. If 5 is convex, then ipQit) € S for all t G [0, 1]. 

The convex hull of a set, V, is the intersection of all the convex sets that contain V. For example 
let J' = {ji,j2, ■ ■ ■ ,jk} denote any fe-element subset of {1, 2, . . . , n + 1}. For ease of notation, we 
assume that ji < j2 < ■ ■ ■ < jk- Then let ej/ = ej^. The vector ej/ is the vector that has 1 
in each of the entries ji and elsewhere. Of course, e^jy = ej, and in this case we use the latter 
notation. By convention, 60 = (0, 0, ... , 0). We consider, V = {eji : J' C {1, 2, . . . , n + 1}}. The 
convex hull of V is the unit cube. Denote the convex hull of any set S, by H{S). 

The construction of the preceeding paragraph indicates the deep connections among the 2"'"'"^ 
vertices of the + l)-cube, the 2"+^ subsets of {1, 2, . . . , n + 1}, and the lattice of inclusions among 
these subsets. Specifically, if J" = J'U {ik+i}, then there is an edge of length 1 in the (n + l)-cube 
that connects the corresponding points ej" and e j' . 

4.2 The {n + 1) -dimensional pyramid 

As a further example, consider any subset J = {ji,j2,---,jk} C {1,2, ...,n}, and let V = 
{e,/u{n+i} ■ J C {1, 2, . . . , n}}. Then V consists of the union of the vertices of the n-cube, 
C"(l) = H{V). The 1 is meant to denote that this cube is the "ceiling" of the (n + l)-cube. 
Now let Py„_|_i denote the convex hull of F U{e0}. The set Py„^i is a "pyramidal" set that consists 
of the set of all lines segments from a point in the "horizontal" n-cubc C"(l) = {(xi,X2, ■ ■ ■ , Xn, 1) G 
: < < 1 for all j = 1, 2, . . . , n} to the origin (0, 0, . . . , 0). The intersection of the pyramid 
Py„+i with the horizontal hyperplane Xn+i = x is an n-cube with the length of each of its edges 
equal to x. Thus the (n + l)-dimensional volume of Py„+i coincides with dx. 

Let us examine in some detail the claim that the intersection of the pyramid Py„+i with 
the hyperplane Xn+i = x is an n-cube with each of the lengths of its edges equal to x. Let 
P = {xi,X2, . . . , Xn, 1) denote a point in the horizontal n-cube, so that for j = 1, 2, . . . , n, we have 
< Xj < 1. The pyramid consists of all lines from the base cube to the vertex Q = (0, 0, . . . ,0). 
Such a line would be 

ipgit) = {tXi,tX2, . . .,tXn,t). 

If t = X, since < x^- < 1, we have < xxj < x, for j = 1,2, ... ,n. That is, each of the early 
coordinates is bounded above by x. 

Now we examine the claim that the (n + l)-dimensional volume of Py„_|_i is Jq x" dx. For each 
X G [0, 1] the height of the curve x" is the n-dimensional volume of the n-cube each of whose edge 
lengths is x. Cavalieri's principle (see any standard calculus text) states that if two figures have the 
same cross sectional measure then the integrals coincide. This idea was used by Archimedes, for 
example, to show that the volume of a sphere could be obtained from the volume of a cylinder and 
the volume of a cone. Here we use this idea to translate the integral x" dx into a computation 
of the volume of the pyramid Py„_(_]^. Upon careful consideration, it seems much more natural to 
consider this rectilinear figure rather than the curved region bounded above by the power function. 
The power function is an expression of n-dimensional volume. The integral adds these volumes to 
get an (n -|- l)-dimensional volume. 
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4.3 An n-dimensional simplex 



Among the pantheon of geometric figures in the plane, the square and the equilateral triangle stand 
out. The 3-dimensional analogue of a square is a cube, and the 3-dimensional analogue of a triangle 
is a tetrahedron. The n-dimensional simplex, or n-simplex, is the higher dimensional analogue of 
the triangle and tetrahedron. It is a regular figure and can be thought of as the controlling figure 
of higher dimensional geometry in the following sense. 

In Euclidean geometry, we learn that any two points determine a line, and any three non- 
collinear points determine a plane. In such a plane, there is a triangle that is the convex hull of 
the three given points. Any four points that are not coplanar determine, via their convex hull, a 
tetrahedron. And that's as far as things go if you limit yourself to exist in the 3-dimensional world. 
If on the other hand, you see that the analogy should carry on, we could consider any five points 
that do not co-exist in a single 3-dimensional hyperplane should form the basic figure, a 4-simplex, 
in 4-space. And so forth. 

In the triangle each pair of points is joined by an edge, and the same is true of the tetrahedron. 
It has six edges, four triangular faces, and four vertices. In an n-simplex every pair of vertices is 
joined. Here is a nice way to draw the 1-dimensional spine (edges) of the n-simplex. Consider the 
set of (n + l)-points arranged evenly spaces on a unit circle. Then draw a line between every pair 
of points. For example, the 4-simplex projects to a regular pentagram (see Fig ??. It is amusing 
and rewarding to draw a circle with a protractor, measure an angle of 360/(n + 1), mark points on 
the circle subtended by multiples of these angles and, draw lines between every pair of points. The 
observant artist will see other lower dimensional simplices emerge from such a process. 



Figure 3: The four dimensional simplex from two standard points of view 
An illustration of an 8-simplex is given in Fig. |^ 

We are interested in a particular regular n-simplex, that is the convex hull of a set of points in 
the (n + l)-cube. The points in question consist of ei = (1, 0, . . . , 0, 0), 62 = (0, 1, ... , 0, 0), through 
e„+i = (0,0,..., 0,1). 

In an n-simplex every pair of vertices are joined by a line segment. The distance from Cj to 
Cj, computed using d{ei,ej) = \J (Axi)^ + (Ax2)^ + • • • + (Axn+i)^) = Similarly, every set of 
3-vertices determines a triangle, every set of 4-vertices determines a tetrahedron, and in general 
every set of /c-vertices determines a (A; — l)-dimensional simplex that is the convex hull of that 
subset. Each of these sub- figures is also regular. For definiteness, we label the n-simplex in which 
we are interested as the convex hull of {ei, 62, . . . , e„+i} 

Sn = H{{ei,e2, . . . ,e„+i}). 
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Figure 4: The set of edges in the 8-dimensionaI simplex 



5 The intersection between Py^+i and Sn 

Before we make the computations of this section, we want to remind you of where we are in this 
process. We have argued that 



where by Vol, we mean the (n + l)-dimensional volume of that figure. We want to compute this 
volume by showing that the pyramid takes up exactly, l/(n + 1) of the {n + l)-cube. To see this, 
we are going to compute the intersection of the pyramid with the simplex, and apply a set of 
symmetries of the simplex to the entire (n + l)-cube. As the symmetry is applied the pyramid will 
intersect another portion of the simplex, and this intersection is an indication of how the pyramid 
and its translations fill the cube. Now we commence. 

Let j denote a fixed integer from 1 to n. Consider the edge £j in the pyramid Py^+i that consists 
of the line segment from to ej + e^+i. In parametric form, £j{t) = t{ej + e„_|_i) for t G [0, 1]. 
Meanwhile, the edge in the simplex that joins ej to en+i is given as Sj{u) = ej + u(e„+i — ej) for 
u G [0, 1]. The intersection between these two lines occurs when ej + u(e„+i — ej) = t{ej + e„_|_i). 
Since {cj : j = 1, 2, . . . , n + 1} is a basis for the vector space M"+^, we have 1 — u = t, and u = t. 
Therefore the intersection point is {ej + en+i)/2. 

Let J = {ji,j2, ■ ■ ■ ,jk} C {1,2, .. . ,n}. Then ejy{„_|_ii, is a vertex of the n-dimensional cube 
C"(l) that is the ceiling of the (n + l)-cube. The line ij{t) = tsju^^^iy = tej + tcn+i is the line on 
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the cone from the origin to that vertex — it is an edge of the pyramid Py„+i . This edge intersects 
the simplex Sn in a A;-dimensional face that is determined by the subset J = {^1,^2, ■ ■ ■ ,3k\- Let 

Tjyj{n+i} = H{{ej^,ej2,. . . , Cj^^, e„+i}). 

Then Tjij^n+i} has the geometric shape of a fc-simplex. For example, if A; = 1, then Ty^ ^+i} is a hne 
segment of length -v/2- If A; = 2, then T^jij.j.n+i} is an equilateral triangle with vertices ej^, ej^ and 
Cn+i- A linear combination, X^j^;^ Uiej^+Un+len+l is an element of Tj provided X^j^j^ Ui + Un+i = 1 
(and Ui>0). 

The intersection between Tj^^^^iy and ij{t) is computed: 

k k 

i=l i=l 

We equate the coefficients of Cj. to get Ui =t for each i = 1, 2, . . . , n + 1. Since Yli=i + '"n+i = !> 
t = + 1) and the intersection point is 

k 

Pju{n+i} = + + + l)en+i = l/{k + 1) X + + l)en+i- 

The reader will observe that this computation used the fact that £ = {ei, 62, ■ ■ ■ , Gn+i} is a basis, 
so that there is a solution to these equations and the solution is unique. 

For J = {ji,j2, . . . , Jfc} C {1, . . . , n+1}, let pj = ^ ^jeJ '^^e intersection between Py„+i, 
and Sn consists of the set 

L>„[n + 1] = H{{pj : n + 1 G J C {1, 2, . . . , n + 1}}). 

The set Dn[n + 1] is an n-dimensional cubiod with its two dimensional faces being kites. It is 
combinatorially equivalent to the ra-cube, C"(l) on the ceiling. More generally, let Dn[i] = H{{pj : 
i e J C {1,2, . . . ,n + 1}}). Each Dn[i] is an n-cuboid, and the union of these for z = 1, . . . , n + 1 
tiles Sn- 



6 Symmetry of the (n + l)-cube 

Recall that the line ^{i,2,...,n}(*) = {t^t-i ■ ■ ■ ^t) is the main diagonal of the (n+ l)-cube. It joins the 
points (0, 0, ... , 0, 0) and (1, 1, ... , 1, 1). We will rotate the cube about the axis 2 i^ such 

a way that ei rotates to 62, 62 rotates to e^, and so forth, e„+i rotates back to ei. The rotation 
which acheives this is given by 



e 



n+l 



/ 

1 

1 
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With @n+i given in this form, for n an even integer, det (6n+i) = —1; for n an odd integer, 
det(G„+i) = 1, and in general {@n+i)^'^^ = Id. In other words 0„+i(xi, X2, . . . , x„, x^+i) = 

Note that G„-|-i is a symmetry of the (n + l)-cube since it preserves the set of vertices. Since 
Sfi passes through the points ei to e^+i, and @n+i preserves the set {ei, . . . , e^+i}, Qn+i is a 
symmetry of Sn. 

We can see that the action of 6n+i on Dn[n + 1] is to move it to Dn[l] as follows. For 
any set J = {ji,j2, ■ ■ ■ ,jk} C {1, 2, . . . , n + 1} and for any integer i, let J + i denote the set 
{ji + h 02 + i, ■ ■ ■ iJk + i} where the elements are reduced modulo n + 1 to lie between 1 and n + 1 
inclusive. Then Qn+i{pj) = PJ+i, and {Qn+iYiPj) = PJ+i- 

7 The link of the vertex 60 and the rotations of the pyramids 

In the {n + l)-cube C„+i, we have a notion of a link of any particular vertex. Here we only 
need to examine the concept for the vertex eg = (0,0, ... ,0). The link i(e0) consists of those 
n-dimensional cubes that do not contain eg. These cubes are the intersections of C„+i with the 
hyperplanes: xi = 1, X2 = 1, . . = 1 and Xn+i = 1. Clearly, eg is not a point in any of these 
planes. Addresses in each of these planes can be determined by the remaining unspecified vectors. 
Thus to find your way around xi = 1, you can use the direction vectors, 62,63, through e^+i, and 
in Xj = 1 you can use ei, through ej_i, and 6j+i through 6„+i. Finally in Xn+i = 1 use 61 through 

6n- 

The rotation @n+i, then rotates these planes among each other. The frame of reference 
(ei, 62, . . . , e„) in Xn+i = 1 rotates to the frame of reference (62, 63, . . . , 6„, e,i+i) in xi = 1. 
Similarly, the frame of reference (62, 63, . . . , 6„, 6„+i) for xi = 1 transforms to the frame 

(63,64, . . . ,6„,6i) for X2 = 1. 

These descriptions of the transformations indicate how the n-cubical base of the pyramid Py^+i 
transforms under the rotations. Since the pyramid is formed by coning this base to eg, this gives 
(n + l)-copies of the pyramid. They fill the (n + l)-cube as we are coning the link ^(eg) to eg. 
Thus Vol(Py„_,_i) = l/(n + 1). In even dimensions, the pyramids under go reflections at each stage 
of the transformation. 

This completes the proof that Jq x"' dx = l/(n + 1) for n > 0. In the next, brief section we 
discuss how the case n = fits into the same framework. In the final section, we indicate how to 
compute the hyper- volumes of some associated n + 1-simplices. 

8 Decomposing Cn+i into simpUces 

Observe that the square [0, 1] x [0, 1] is cut into two triangles via the pyramid [sic triangle] Py2 and 
its reflection. We can use that decomposition, in the pyramid Py3 to first decompose the base into 
two triangles, and subsequently decompose the pyramid into two tetrahedra. One is the reflection 
of the other; they are obtained by "coning off" the diagonal of the ceiling to the point eg = (0, 0, 0). 
In this way we can cut the cube into six congruent pieces, each of which is a tetrahedron. The base 
of the tetrahedron has area 1/2. the height of each is 1. The volume, consequently is 1/6. 
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Figure 5: The link of (0, 0, 0, 1) consists of four cubes of dimension 3 

The process continues by induction. The (n + l)-cube can be decomposed into a total of (n + 1)! 
congruent (n + l)-simplices, each of which has hyper- volume l/(n + 1)! The construction follows 
by decomposing the base cube of the pyramid Py^+i into n! simplices, coning within the pyramid, 
and rotating using 6. 

This decomposition is interesting for some historical reasons. It is impossible to cut a 3-cube of 
volume 1 into finitely many pieces, and reassemble the result into a regular tetrahedron of volume 
1. Such a process is known as a "scissors congruence." The impossibility was demonstrated by Max 
Dehn who used an invariant defined via the dihedral angles of the cube and tetrahedra. Scissors 
congruences give important invariants of 3-dimensional solids. We remark that AC studies these. 

JSC uses decompositions of higher dimensional solids to identify coincidences among certain 
cohomology classes. These cohomology theories are used to develop invariants of higher dimensional 
knots, and the cohomology theories are used to describe and analyze algebraic structures. 

Therefore, this paper illustrates the value in "thinking deeply about simple things." 
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9 Analytic description 



This tesellation of the (n+ l)-cube has a very pretty analytic description. For each i = 1, . . . , (n+ 1) 
let 

Pyn+l = {(a^l' ■ ■ -^Xn+l) e Cn+l\ Xj < Xj Vj = 1, . . . , (n + ^ i} 

Tha plane Xi = x intersects Py^+i in the n-cube with the length of each of its edges equal to x. 
Hence Vol(Pyl+i) = x^dx. Also G„+i(Pyj,+i) = Pyj+\ for alH = 1, . . . ,n and @n+i{^jlX\) = 
P/n+i hence {Py^+i| i = I,. . . , (n+1)} tesellate the (n+l)-cube showing dx = l/(n+l) 

for n > 

10 Mathematica code 

The code below was writing with anearlier version of the paper in mind. In the earlier version the 
cube of size x^ was located in the plane Xn+i = 1 — a; so that the pyramid Py„+i had its base on 
the floor of the (n + l)-cube. That turned out to make the rotations of the cube appear to be very 
complicated. On the other hand, it is not a great idea to modify computer code (even inelegant 
computer code) that runs. So there are inconguencies between the test and the program. 



vera[prop_] := If [TrueQ [prop] , 1, 0]; 
val [t_ , a_ , b_] : = vera [a <= t] vera [t < b] ; 
vem[t_, a_, b_] := vera [a <= t] vera[t <= b] ; 
lyle[n_, j_, t_] := Binomial [n, 2] t - j +1; 

booze[n_] := Flatten [Table [{i, j>, {i, 1, n - 1>, {j, i + 1, n}] , 1]; 

wine[j_, n_, t_] := ReplacePart [ 

ReplacePart [ 
ReplacePart [ 

ReplacePart [IdentityMatrix[n] , 

Cos[ 2Pi lyle[n, j, t] ], {booze [n] [ [j] ][ [1] ] , 
booze[n][[j] ][[!]]>], 

Cos[2Pi lyle[n, j, t]], {booze [n] [ [j] ] [ [2] ] , 
booze [n] [[j]] [[2]]}] , -Sin[ 

2Pi lyleCn, j, t]], {booze [n] [ [j] ] [ [2] ] , booze [n] [ [j] ][ [1] ]>] , 
Sin[ 2Pi lyle[n, j, t]], {booze [n] [ [j] ][ [1] ] , booze [n] [ [j] ] [ [2] ] }] ; 
tire[t_, n_] := 

Suin[val[t, (j - 1) /Binomial [n, 2], j/Binomial [n, 2]]wine[j, n, t] , { j , i, 
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Binomial [n, 2] - 1}] + 
vein[t, (Binomial [n, 2] - 1) /Binomial [n, 2], 1] wine [Binomial [n, 2], n, 
t]; 

(*This cell establishes a set of rotations in each of n choose 2 directions, 
one performed after the other. 

The rotations can be appied to any higher dimensional figure *) 



ramona[m_] :=Table [1/Sqrt [m] ,{i,0,m-l}] 
(* ramona is used to give a 3- 

dimensional projection. It may not be the best choice *) 

Clear [eh] (* eh is Canadian, 

eh? Now you have to be careful how you choose to execute the 
next few cells. *) 

eh[4]=Transpose[{{0,l},{l,l>/Sqrt[2] , {1 ,0> , {1 , -1}/Sqrt [2] }] (* 
this value of eh is good to project the 4- 
cube into the standard octahedral projection. 
It is not much good for any thing else. *) 

eh[m_] :={Table[Re[E'-(2 Pi j I/m)] ,{j ,0,m-l}] , 

Table[Im[E-(2 Pi j I/m)] ,{j ,0,m-l}] } 
(* This version of eh gives a 2dimensional projection to the plane where the \ 
coordinate axes go to the roots of unity. 

These projections tends to give the most symmetric pictures. *) 

eh[m_] :={Table[Re[E~(2 Pi j I/m)] ,{j ,0,m-l}] , 

Table[Im[E"(2 Pi j I/m)] ,{j ,0,m-l}] ,ramona[m] } 
(* This version of eh projects to 3 dimension, 
it includes some strange distortions. *) 



eh[m_] :={IdentityMatrix[m] [[1]] ,IdentityMatrix[m] [[2]] , 
IdentityMatrix[m] [[3]]} (* 
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The most obvious choice for eh. 

You can use this to gain more intution since various details are \ 
squished out . *) 

In [8] : = 

eh[m_] :={ldentityMatrix[m] [[1]] ,IdentityMatrix[m] [[2]] , 
Identity-Matrix [m] [ [3] ] +IdentityMatrix [m] [ [m] ] > (* 
Another obvious choice for eh. Try running with this value. 
This is my favorite view for the 4- 
cube. You really can see the cuboid in the simplex here 
even though other details are squashed out.*) 



(* CAUTION. ONLY EXECUTE ONE VALUE OF eh. 

ALSO YOUR CHOICE OF eh DETERMINES WHETHER 
YOU RUN moonbeam or \ 
moonbeam2d! *) 

In [9] := 

topvert[n_] :=Insert [Table [0,{i,l,n-l}] ,l,n] ; 

missing [j_,m_] :=Join[ Take [Array [x,m] ,j-l] , Take [Array [x , m] ,j-m]] ; 
nerd [j_ ,m_ , s_] :=Insert [missing [j ,m] ,s, j] ; 
felix[j_,m_,t_] := 

Line [{Dot [eh [m] ,tire [t ,m] ] .nerd [j ,m, 0] , 
Dot [eh [m] ,tire [t ,m] ] .nerd[j ,m, 1] }] ; 
leonard[j_,m_,s_] :={Hue[.0] , 

Line [{Dot [eh [m] ,tire[t,m]] .nerd[j ,m,s] , 
Dot [eh[m] ,tire [t ,m] ] . topvert [m] }] } ; 
git [m_] : =Table [IntegerDigits [j , 2 ,m] , {j ,0 , 2"m-l}] ; 
shirt [m_] : =Table [IntegerDigits [2j ,2,m] ,{j , 0,2~m-l}] ; 
francos [j_ ,m_ ,k_ ,git_] : = 

Table [missing [j ,m] [ [i] ] ->git [m-1] [[k,i]] ,{i,l. 
Length [missing [ j ,m] ] }] ; 
hyper [j_,m_,t_] := 

Table [f elix [j ,m, t] / . f ranees [j ,m,k,git] , {k , 1 , Length [git [m-1] ] }] ; 
hypo[j_,m_,t_] :={Hue [j/m'"2] , 

Table [f elix [j ,m,t] / .f ranees [j ,m,k, shirt] , {k, 1 .Length [shirt [m-1] ] /2}] } ; 
(*This cell provides the main portion of the hyper cube, topvert is e_{n+l}. 
Felix and leonard are lines. 

Most anything with a feline sounding name is a line: line-> 
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lion->cat. git stands for digit, 
shirt is the sanitized version of a rhyming word, 
f ranees is a substitution rule: rule-> 
mule-> 

f ranees the talking mule. It is misspelled on purpose, 
hypo and hyper are the corresponding cubes . 
The edges of hypo are specially colored. *) 



In [19] : = 

parrot [ j _ , k_ , m_ , 

t_] :={Surf aceColor [ 

RGBColor [1-Sqrt [ ( j ~2+k~2) ] / (2m) ~2 , , 1- j /ni'2] ] (*Hue [ ( j+k) /4m] *) , 
Polygon[{Dot [eh[m] ,tire [t,m]] . Array [x,m] /. {x [j] ->0 ,x [k] ->0>, 
Dot[eh[m] ,tire[t,m]] . Array [x,m] /. {x [j] ->1 ,x [k] ->0}, 
Dot [eh[m] ,tire[t,m]] .Array[x,m]/.{x[j]->l,x[k]->l}, 
Dot [eh[m] ,tire[t,m]] . Array [x,m] / .{x[j] ->0,x [k] ->1}}] }; 
twoholes [j_ ,k_ ,m_] : = 

Join[Join[ Take [Array [x , m] ,j-l] , Take [Array [x , m] ,{j+l,k-l}]] , 
Take [Array [x , m] , k-m] ] ; 
f ools [j_,k_,m_] :=Table [twoholes [j ,k,m] [ [ell] ] ->0 , {ell , 1 ,m-2>] ; 
rooster [j_ ,k_ ,m_ ,t_] :=parrot[j , k,m,t]/. fools [j ,k,m] ; 
hone [j_ ,m_ ,t_] := 

Table [leonard[j ,m,0] / . f ranees [j ,m,k,git] ,{k, 1 , Length [git [m-1] ] }] 
(* Things avian are surfaces: It's a bird its a plane, hone is the cone. 
I use the Dr. 

Suess approach to programming to avoid Mathematica thinking that \ 
uncapitalized words are misspelling of its commands . 

fools are a set of substitution rules.*) 

In [24] : = 



jag[u_,v_,m_,t_] :={Hue[0.5] , 

Line [{Dot [eh [m] , tire [t ,m] ] . (IdentityMatrix [m] [ [u] ] +topvert [m] ) , 
Dot[eh[m] ,tire[t,m]] . (IdentityMatrix [m] [[v]]+topvert [m] )}]}; 
uar [u_ ,m_ ,t_] :={Hue [0.5], 

Line [{Dot [eh [m] ,tire[t,m]] . (IdentityMatrix [m] [ [u] ] +topvert [m] ) , 
Dot [eh [m] ,tire[t,m]] . (IdentityMatrix [m] [ [m] ] -topvert [m] ) }] } ; 
pimplex[m_,t_] :={AbsoluteThickness [3] , 

Join[Table [jag[u,v,m,t] ,{u,l,m},{v,u+l,m-l}] , 
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Table [uar [u,m,t] ,{u,l,iii-l}]]} 
(* jag and uar are sets of lines, pimplex is the simplex.*) 



In [26] := 

powder [iii_] : =Table [Reverse [IntegerDigits [k,2,iii] ] , {k, 2~ (m-l)+l , 2~m-l}] 
(* powder was my approach to finding the power set. 
This portion of the program is 
very unelegant . *) 



In [27] := 

bims [k_,m_] : = 

Sum [powder [m] [ [k, ell] ], {ell , 1 ,m}] (* buns, sums, 
buns are made of flour a type of powder *) 

In [28] := 

mist [m_] : =Do [spot [ j ,m] ={} ; 

Do [If [Equal [buns [k,m] , j] , AppendTo [spot [j ,m] .powder [m] [ [k] ] ] ] ,{k, 1 , 
Length [powder [m] ] }] , {j , 2 ,m}] ; 
(* mist is a given list. This is a very weird peice of programming. 

See the next comment *) 



In [29] : = 

Table[mist[k] ,{k,3,7}] (* 

OK YOU HAVE TO RUN THIS FOR THE LAST NUMBER AT THE DIMENSION YOU WANT \ 
TO GRAPH. THE PROGRAM WORKS IN THEORY FOR ALL N, 

BUT ON MY MACHINES I NEVER GOT PAST 7. I NEED A NEW COMPUTER! IGNORE THE \ 
NULL OUTPUT.*) 

Out [29] = 

{Null , Null , Null , Null , Null} 
In [30] := 
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Clear [freckle] 
In [31] : = 

freckle[k_,iii_,ind_] := 1/k spot [k,m] [ [ind] ] + (k-2)/k IdentityMatrix [m] [ [m] ] 
(♦freckle [k.ni, ind] is a barycenter of the k-1 dim'l simplex in m- 
space . The index= 
ind is a numbering thereof. 2\ [LessEqual]k\ [LessEqual]m, 
and ind\ [LessEqual] 
Binomial [m-1 , k-1] . 

The freckle is the spot where the simplex intersects the cone *) 

In [32] := 

Clear [victor] 
In [33] := 

victor [k_,ell_,m_,ind_,dex_] : =spot [ell ,m] [ [dex] ] -spot [k,m] [[ind]] 

(* What's the vector victor? What's the clearance, Clarence? Airplane *) 



In [34] := 

Clear [guess] ; Clear All [dorothy] 

In [35] := 

guess [m_] := 

Do [dorothy [k,m] ={]■ ; 

Do [If [victor [k,k+l ,m, ind, dex] \ [Equal] Abs [victor [ k,k+l ,m, ind,dex] ] , 

AppendTo [dorothy [k,m] , {ind, dex}] ] , {ind, 1 , Binomial [m-1 , k-1] } ,{dex, 1 , 
Binomial [m-l,k]}] ,{k,2,m-l} ] 
(* dorothy is supposed to be close to the cowardly lion. 

dorothy gives the set of lines of intersections between the faces \ 
of the cone and the simplex. *) 

In [36] : = 

Table [guess [m] ,{m, 2, 7}] (* 

I guess I got this part right. It looks good in low dimensions. 
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THE LST NUMBER SHOUID BE EQUAL TO OR EXCEED THE DIMENSION YOU WANT TO \ 
GRAPH! BAD PROGRAMMING *) 

Out [36] = 

{Null , Null , Null , Null , Null , Null} 

In [37] := 
Clear [cello] 

In [38] : = 

cello [k_,ell_,m_,ind_,dex_,t_] :={AbsoluteThickness [2] , 
Hue [0.7] (*Hue[l-(k+ell)/m-3]*) , 
Line [{Dot [eh[iii] ,tire [t.m]] . freckle [k,m,ind] , 

Dot[eh[m] ,tire[t,ni]] . freckle [ell, m,dex] }] }; 
(* I think I wanted violins here, 
but I used a cello instead. I can't remember why. 
You can change the Hue if you like . 
I liked .7 since it was a nice blue.*) 

In [39] := 



vegetable [k_,ell_,m_,t_] := 
Partition [ 

Flatten [Table[ 

cello [k, ell ,m,dorothy [k,m] [ [j , 1] ] ,dorothy [k,m] [ [j ,2] ] ,t] ,{j , 1 , 
Length [dor othy[k,m]]>]] ,3] (* 
vegetables such as carrot sticks. These are the sticks in the simplex*) 

In [40] := 

salad [m_,t_] := 

Table[vegetable[k,k+l,m,t] ,{k,2,m-l]-] (* a bimch of vegetables *) 
In [41] : = 

polka [k_,m_,ind_,t_] :={PointSize[0.01] ,Hue[0.99] , 

Point [Dot [eh [m] ,tire [t ,m] ] . freckle [k,m, ind] ] }; 
(* a bunch of dots *) 

In [42] : = 

dance [m_,t_] := 
Partition [ 
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Flatten [ Table [polka [k,m,ind,t] , {k, 2, m}, {ind, 1 .Length [spot [k,m] ]}] ] ,3] 
(* such as a polka *) 

In [43] : = 

ohmy[m_,t_] := {PointSize [0 . 02] .Point [Dot [eh [m] ,tire[t,m] .Table [0,{j ,l,m}]]]} 
(* the origin of the cube indicated in a big black dot. *) 



(* RUN THIS IF eh is 2 by m *) 
moonbeaiii2d [m_ , dt_ , a_ , b_] : = 
Table [Show [ 

Graphics [{Table [{hyper [j ,m,t] } , {j , 1 ,m}] .{AbsoluteThickness [2] , 
Table[{hypo [j ,iii,t] >, {j , 1 ,iii-l>] >, 
hone [nijiiijt] , 

salad [m,t] , dance [m,t] ,pimplex[m,t] ,ohiny[m,t] }] , 
PlotLabel-> 

StyleForm[IntegerPart [t Binomial [m,2] /dt] , "Section" , 
FontColor->GrayLevel [ . 1] ] , AspectRatio->l , Axes->False, 
IinageSize-> {(40*6) , (40*6) } (* ,BoxRatios->{l , 1 , 1}, 
Boxed->False*)] ,{t ,a,b,dt/Binomial [in,2] }] 
(* m is the dimension (n+1) in the paper. 

dt is an increment of the movies a and b should be between and \ 
1 . you can control these if you want one rotation *) 

moonbeam2d [4,0.0125,0,1] 

(* THE LINE THAT GIVES OUTPUT IN CASE eh is 2 by m *) 
In [44] := 

(* RUN THIS IF eh is 3 by m *) 
moonbeam [m_ , dt_ , a_ , b_] : = 
Table [Show [ 

Graphics3D [{Table [{hyper [j ,m,t] } , {j , 1 ,m}] , {AbsoluteThickness [3] , 
Table [{hypo [j ,m, t] } , {j , 1 ,m-l}] } , {AbsoluteThickness [2] , 
hone [m,m,t] }, Table [rooster [j ,k,m, t] , {j , 1 ,m-l} ,{k, j + 1 ,m-l}] , 
salad [m,t] , dance [m,t] ,pimplex [m,t] ,ohmy[m,t]}] , 

PlotLabel-> 

StyleForm[IntegerPart [t Binomial [m,2] /dt] , "Section" , 
FontColor->GrayLevel [ . 1]] ,AspectRatio->l,Axes->False, 
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ImageSize-> {(72*6) , (72*6) > ,BoxRatios->{l , 1 , 1} , 

Boxed->False] , {t , a, b,dt /Binomial [m, 2] }] 
(* m is the dimension (n+1) in the paper. 

dt is an increment of the movies a and b should be between and \ 
1. you can control these if you want one rotation *) 

moonbeam [4, .0125,0,1] (* THE LINE THAT GIVES OUTPUT IN CASE eh is 3 by m *) 

11 A gallery of imagery 

Our inclination in this section is to provide a collection of images that were produced by minor 
variants on the mathematica code provided. The variations occur in the choices used for the 
projection matrix 



Not every choice that we used are listed in the program code. However, we include some minor 
comments. 

In the first illustration, the figure resembles the drawn image. The reason is that in both cases 
the octohedral view point was chosen. The rough sketches of the illustration were done with the 
computer off, and the Mathematica illustration was used to check that coordinates were correctly 
indicated. The cyan in the drawing was chosen to coincide with the color choice in Mathematica. 
The next two illustrations are from the same movie. 



eh. 
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The next figure helps the reader see that the intersection of the pyramid with the tetrahdron 
is a cuboid. 

180 




This ihustration was chosen using the squished view of the four simplex. The tetrahedron is 
squashed into a triangle, but the cone to the floor is quite visible. 

430 
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Here is one nice view of the 5-cube with cone simplex and intersection indicated. 



137 




This view of the six dimensional cube uses a different value for 
eh. 

On our machines, we can run the corresponding illustration with a large number of stills per 
movie, but then we cannot animate the result. 






In principle, the program code listed below runs to illustrate any n-cube. In practice it gets 
bogged down at about n = 6. The top of the code can be used to illustrate any higher dimensional 
object. The object is rotated full circle in each of the standard (2) planes in n-space. Within the 
hyper-cubes, only certain 2-dimensional faces are shown again of them. 
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Here is a fake stereopsis that was created by taking two successive views of the 4-cube and 
setting them side-by-side. Try crossing your eyes to see it in 3D. 



387 
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